Large Contributions of Negative Energy States to Forbidden Magnetic-Dipole 
Transition Amplitudes in Alkali-Metal Atoms 
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The influence of negative-energy states (NES) on forbid- 
den magnetic-dipole nsi/ 2 — (n + l)si/2 transitions in alkali- 
metal atoms is investigated. We find that the NES contri- 
butions are significant in almost all cases and for rubidium 
reduce the transition rate by a factor of 8. We tabulate 
magnetic-dipole (Mi) transition amplitudes for the alkalis. 
Our Mi value for cesium, where accurate measurements are 
available, differs from experiment by 16%. We briefly discuss 
the feasibility of an experimental test of NES effects. 

PACS number(s): 31.30. Jv, 31.15.Md, 32.70.Cs 



It is well known (see, for example, the discussion 
in Ref. [Q) that the Dirac-Coulomb Hamiltonian has 
no bound state eigenfunctions in the presence of the 
electron-electron interaction. The no-pair Hamiltonian 
H n .p. derived from quantum electrodynamics has been 
proposed by Brown and Ravenhall Q for use in relativis- 
tic atomic calculations. Although H n . p . leads to very ac- 
curate energies, the omitted effects of electron-positron 
pairs can be significant for the forbidden magnetic-dipole 
(Mi) transition amplitudes. The first discussion of pair 
corrections to the Mi decay rate for the 2 3 S\ state 
in helium was given by Feinberg and Sucher || in 
1971. Later, when new lifetime measurements of the 
2 3 Si metastable state of heliumlike ions became avail- 
able, more accurate calculations of Mi were performed 
by several theoretical groups. Lindroth and Solomon- 
son ||] numerically demonstrated the detailed cancella- 
tion of one-pair diagrams in the case of heliumlike ar- 
gon. The decay rates of the same transitions for heli- 
umlike ions with Z = 2 — 100 were calculated within 
relativistic configuration-interaction approach by John- 
son, Plante, and Sapirstein ||. They treated contribu- 
tions of NES using second-order many-body perturbation 
theory (MBPT). Indelicato || considered the effects of 
NES for 2 3 S\ Mi decay in the multiconfiguration Dirac- 
Fock (MCDF) approach. Recently we studied NES con- 
tributions to transition amplitudes in more detail @. 
In the Pauli approximation, we derived an effective one- 
pair operator that explicitly reveals cancellation between 
Coulomb and Breit two-body diagrams and, by using it, 
found a transition without such cancellation: the neu- 
tral helium 2 3 Si — 3 3 Si transition has a very large NES 
contribution which reduces the Mi rate by a factor of 
2.9. 

There have been no other calculations treating NES 
contributions to Mi transitions systematically in multi- 



electron systems, except for He- and Be-like transi- 
tions ||. In this letter, we report second-order MBPT 
calculations of magnetic-dipole nsi/ 2 — (n + l)si/ 2 tran- 
sitions in the alkalis including the analysis of one-pair 
effects. We have discovered an unusually large NES con- 
tribution, which reduces the transition rate by a factor of 
8 in Rb, and propose measurements to test NES effects. 
The results of our calculations for Cs are compared to a 
previous theoretical determination Jj| and to experimen- 
tal values |l(],[ll|]. For the other alkalis, no measurements 
exist. 

We can argue that forbidden magnetic-dipole ampli- 
tudes are the most sensitive among electro-magnetic 
transition amplitudes to the accuracy of the relativis- 
tic description of an atomic system. As we will demon- 
strate, several factors contribute to the result: correlation 
effects, spin-orbit interaction, Breit interaction, retarda- 
tion effects, and, finally, the negative-energy contribu- 
tions. The relativistic retarded magnetic-dipole matrix 
element can be represented as (5| (atomic units are used 
throughout the letter) 



(z||Mi||j) = C (-Ki\\Ci\\Kj) (Ki + Kj) X 
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(GiFj+FiGj) dr. 



(1) 



Here G and F are the radial parts of large and small com- 
ponents of atomic orbitals, k is the photon wavenumber, 
and k = (j + i) (— iy+ l + 1 / 2 . In the long-wavelength 
limit and Pauli approximation this relativistic expression 
reduces to a conventional non-relativistic operator 



Mi = L + 2S . 



(2) 



Even if the general angular selection rules for the Mi op- 
erator are satisfied, this matrix element vanishes when 
the radial wavefunctions are orthogonal, i.e. if Ki = Kj, 
but rii 7^ riy The Einstein ^4-coefficient for the Mi tran- 
sition \I) — ► \F) is expressed in terms of the reduced 
matrix element as 



A Ml — 



k 3 \{F\\Mi\\I) \ 
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We start our consideration by utilizing second-order 
MBPT built on the "frozen-core" Dirac-Hartree-Fock 
(DHF) potential. This approximation includes both lead- 
ing correlation and NES effects. We consider matrix el- 
ements of the magnetic-dipole operator z between two 
valence states v and w. For the purposes of this pa- 
per, the valence state v represents the ground state or- 
bital nsi/2 and the state w represents the first excited s 
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state (n + l)si/ 2 - The first order value is given by the 



matrix element taken between DHF orbitals z v 



The 



second-order correction including both Coulomb (gijki) 
and Breit (bijki,bij) interactions is given by 



z& = V 



E 



na 



(fju 



e a + £v 



(3) 



Here antisymmetrized matrix elements are defined as 
gijki = 9i 3 ki - 9i 3 ik- The one-body matrix element of 
the Breit interaction is 6y = ^2 a bi a ja- In the above ex- 
pressions, index a runs over core orbitals and n extends 
over virtual orbitals. It is important to note that virtual 
states include both excited positive and negative-energy 
orbitals. Corresponding Feynman diagrams are given in 
Fig. [I]. We emphasize that in the no-pair approximation 
the summation over virtual states would be limited only 
to the positive-energy states, i.e. include only the dia- 
grams in the upper panel of Fig. [l]. The inclusion of the 
Breit interaction in our analysis is crucial, because its 
contribution to the matrix element is of the same order 
as for the Coulomb interaction. We use a static limit of 
the Breit interaction, since the energy of the transverse 
photon is determined by the energy difference of the real 
electrons, not by that of the virtual electron. The nu- 
merical calculations were performed with a relativistic 
B-spline basis set representation fl2]| obtained in a cavity 
of radius 40 atomic units and included 40 positive- and 
40 negative-energy wavefunctions for each partial wave. 

The results of our calculations are presented in Table |. 
Note that the first-order Dirac-Hartree values dominate 
for lighter atoms and become less significant for cesium 
and francium. This is due to larger second-order no-pair 
contributions. The values of NES contributions (in the 
third row) appear to be roughly proportional to Z. In 
the case of cesium they constitute only a small fraction 
4% which is even smaller for francium 0.6% due to large 
no-pair second-order contributions. However, the NES 
fraction reaches 19% in potassium. The rubidium case 
is the most surprising: there is cancellation of the two 
no-pair terms and consequently a strong dependence of 
the total value on the negative-energy corrections. Al- 
though such cancellation in second order may be coinci- 
dental, and more accurate calculations may be necessary, 
we conclude that a measurement of the Mi transition in 
rubidium could provide an excellent test of the NES con- 
tributions. 

The large relative contribution of NES for forbidden 
magnetic-dipole transitions is caused by several factors. 
First, due to unique properties of the Mi operator (||), 
the no-pair amplitude is severely suppressed (by a fac- 
tor of a 2 in the lowest order). Second, the magnetic- 
dipole operator (Mi)jj in Eq. ([!]) contains an integral 
of large and small components of Dirac wavefunctions. 
For positive-energy states the small component is signif- 



icantly weaker than the large component (by a factor of 
aZ for hydrogen- like ions). For NES, the situation is 
reversed, and the small component is much larger. In 
addition, the Pauli approximation expression (|2|) with 
its particular (5-function-likc properties is no longer valid 
and one obtains much larger values for Mi matrix el- 
ements between negative and positive states, than from 
positive-positive matrix elements. Finally, the energy de- 
nominators of order 2mc 2 bring the NES contributions 
to the same level as the contribution from the "regular" 
positive-energy states. As seen from Table 0, NES con- 
tributions from the Breit interaction are comparable to 
those from the Coulomb diagrams because the Breit op- 
erator mixes large and small components. 

We note that for Rb, Cs and Fr, correlation effects 
are very important leading to contributions larger than 
the lowest-order DHF values. The mechanism has been 
discussed by Dzuba et al. ||. In the Pauli approxima- 
tion, the Mi matrix element is proportional to the inte- 
gral of the product of the large components between the 
states involved. In the first-order forbidden transitions 
nsi/2 — (?i+l)si/2 between states with different principal 
quantum numbers, the radial wavefunctions are orthog- 
onal and the Mi rate is zero. Although it is not zero 
beyond Pauli approximation, it is strongly suppressed. 
The situation is quite different for P1/2 and p 3 / 2 matrix 
elements which are non-zero due to overlapping radial 
wavefunctions caused by the spin-orbit interaction. As 
a result, the second order contributions dominate due to 
such matrix elements connecting core and exciting states. 
This correlation effect becomes overwhelming for heavier 
elements where spin-orbit coupling is important. 

In Table Ej , we compare our cesium results for the 
magnetic-dipole reduced matrix clement with calcula- 
tions of Dzuba et al. Q and with measurements from 
several experimental groups. The transition amplitude 
used in || is related to the reduced matrix element ex- 
pressed in atomic units as 

(Mi)ampi. = -7= {n w Sii 2 \\Mi\\n v Sii 2 ) x | — |. 

VD c 

The experimental entries for the Mi matrix element in 
Table [II were obtained from measurements of M^ /Mi 
and a semiempirical value [p. if of the off-diagonal hyper- 
fine mixing amplitude M^ = 0.8094(20) x 10" 5 |^b/c|. 
The result of our work, despite approximate treatment 
of correlation effects, is in reasonable agreement (16%) 
with the experimental results. Since the negative-energy 
effects are marginally smaller than these deviations, it is 
not possible to draw definitive conclusions about NES ef- 
fects from available experiments in cesium. The second- 
order expression (|3|) is a leading term of the random- 
phase approximation (RPA). The calculations of Dzuba 
et al. ^ implicitly included the effect of negative-energy 
states due to the reduction of RPA-like diagrams to the 
form of a differential equation. However, their analysis 
did not take into account the Breit interaction. Such an 
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approach misses an important negative-energy contribu- 
tion. Indeed, we demonstrate in Table || that the NES 
contribution from the Breit interaction is much larger 
than that arising from the Coulomb interaction. 

The theoretical calculations of the Mi transition am- 
plitudes in the alkalis clearly demonstrate the important 
role of negative-energy states. We now discuss experi- 
mental possibilities to test these contributions. We com- 
pare the NES fractional contributions, defined as the ra- 
tio of NES to no-pair contributions, in different alkali- 
metal atoms in Fig. ||. In the light alkalis (Li, Na, K) the 
effect is proportional to Z and is maximal for K (19%). 
For heavy atoms such as Cs and Fr, it is small because 
of large no-pair contributions. Rubidium, in the mid- 
dle, has a very large relative effect (65%) and is the most 
promising. If measurements in the other alkalis reach the 
precision achieved in Cs, then all alkalis except Fr will be 
good candidates for testing NES contributions. The ac- 
curacy of the calculations, on the other hand, can impose 
even more severe restrictions than experiment. The ac- 
curacy of our calculations, as seen in the deviation from 
the experiment for Cs, is about 16%; it is expected to be 
better for lighter elements. More accurate (1%) no-pair 
calculations are possible, for example, in the relativistic 
single-double approximation p3| . For Li, precise no-pair 
configuration-interaction calculations [Q are also possi- 
ble with an accuracy much better than 1%. 

In conclusion, wc have presented the results of the 
second-order MBPT calculations for the forbidden Mi 
transitions in the alkalis. Comparisons with experimen- 
tal and theoretical data for cesium have been made. We 
have found very large negative-energy state contributions 
to the Mi transition amplitudes in the alkalis. The NES 
amplitude is dominant in the case of rubidium, which 
could provide the best experimental test of negative- 
energy contributions in atomic structure. 

We thank Savely Karshenboim for useful discus- 
sions and Eugene Livingston for his comments on the 
manuscript. The work of AD and WRJ was supported 
in part by NSF Grant No. PHY 99-70666. 

TABLE I. Contributions to reduced matrix elements of 
the Mi operator in atomic units multiplied by a factor of 10 5 . 
Row 1, lowest-order DHF value; row 2, second-order no-pair 
contribution; row 3, negative-energy state contributions in 
second order; and row 4, total value of Mi matrix element. 





Li 


Na 


K 


Rb 


Cs 


Fr 


z 


3 


11 


19 


37 


55 


87 


I 


0.91 


1.16 


1.15 


1.38 


1.51 


2.09 


II, no-pair 


0.12 


0.03 


-0.08 


-1.86 


-10.69 


-116 


II, NES 


0.02 


0.13 


0.20 


0.31 


0.40 


0.64 


Total 


1.05 


1.06 


1.27 


-0.17 


-8.78 


-113 



TABLE II. Breakdown of negative-energy state contribu- 
tions to the reduced Mi matrix element in atomic units, mul- 
tiplied by a factor 10 5 . Column "Coulomb" represents contri- 
butions from diagram (d), column "Breit two-body" — from 
diagram (eh and column "Breit one-body" — from diagram 
(f ) of Fig. 





Coulomb 


Breit 

two-body 


one-body 


Total 


Li 


-0.015 


0.067 


-0.029 


0.024 


Na 


-0.020 


0.106 


0.047 


0.133 


K 


-0.022 


0.112 


0.106 


0.197 


Rb 


-0.025 


0.154 


0.174 


0.303 


Cs 


-0.026 


0.183 


0.239 


0.395 


Fr 


-0.035 


0.221 


0.450 


0.636 



TABLE III. Comparison of theoretical and experimental 
values for Cs 6s — 7s reduced magnetic-dipole matrix element 
in atomic units, multiplied by a factor 10 . The experimental 
errors are given in parentheses. 



Reference (6gjjMjJj7g) x 10 5 

Theory 

This work -8.78 

Dzuba et al. §, 1985 -13.7 

Experiment 

Bennett and WiemanlQ, 1999 -10.40(0.03) 

Bouchiat and Guena 1988 a -10.5 (0.1) 



a The average of previous experimental results corrected for 
the electric-quadrupole contribution. 
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(a) 



(b) 



(c) 



(f) 



(d) (e) 
FIG. 1. Principal Feynman-like time-ordered diagrams 
contributing to Mi amplitude in the second-order. The wavy 
lines represent photons and solid lines represent atomic elec- 
trons. Double vertical solid line designates inactive (observ- 
ing) electrons. Diagrams (a) and (d) are due to Coulomb 
interaction, (b) and (e) due to two-body static Breit interac- 
tion, and (c) and (f ) due to one-body static Breit interaction. 
Upper panel of diagrams represents no-pair contributions, and 
the lower panel — contributions from negative-energy states. 
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FIG. 2. The relative contributions to the magnetic dipole 
(Mi) matrix element ns- (n + l)s in alkali atoms: the ratio of 
the NES contributions (row 3 of Table to the total no-pair 
contributions (sum of rows 1 and 2 of Table III ) . 



4 



